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1.  GOALS 


The  objective  of  this  contract  is  to  obtain  ionospheric  measurements  from  a  wide  range  of 
geographic  locations  and  to  utilize  the  resulting  databases  to  validate  the  theoretical  ionospheric 
models  that  are  the  basis  of  the  Parameterized  Ionospheric  Specification  Model  (PRISM)  and  the 
Ionospheric  Forecast  Model  (IFM). 

In  the  past  year,  we  have  supported  these  goals  with  the  following  activities: 

•  tested  the  IFM  model. 

•  initiated  the  use  and  production  of  the  Coupled  Ionosphere  Thermosphere 
Electrodynamic  Forecast  Model  (CITEFM). 

•  developed  an  alternate  derivation  of  the  Alfven-Falthammer  formula  for  the  parallel  E- 
field  in  upward  auroral-current  regions. 

•  contributed  to  validation  studies  of  PRISM. 

•  analyzed  the  ionospheric  response  to  the  Great  Magnetic  Storm  of  July  2000. 

•  initiated  validations  of  the  Global  Positioning  System  (GPS)  position  error  maps. 

•  initiated  development  of  a  PRISM  pre-processor  that  will  enable  data  from  various 
sensors  to  be  integrated  into  the  existing  technology  of  PRISM. 

These  activities  resulted  in  three  presentations  at  various  ionospheric  meetings  and  one  published 
journal  article.  Details  of  these  activities,  presentations  and  the  journal  article  are  provided  in 
this  report. 


2.  PROGRESS 

2.1  The  Ionospheric  Forecast  Model  (IFM) 

The  IFM  model  is  a  collection  of  programs  that  provide  a  first-principles,  global  representation 
of  the  ionosphere  given  a  set  of  geophysical  conditions.  In  the  current  year,  the  IFM  model  was 
utilized  to  generate  ionospheric  characteristics  for  three  levels  of  magnetic  activity  and  three 
levels  of  solar  activity.  The  Total  Electron  Content  (TEC)  and  maximum  density  of  the  F2 
region  (NmF2)  output  from  the  nine  IFM  runs  were  graphically  animated  and  displayed  on  a  web 
site  maintained  by  Boston  College.  Producing  these  animations  facilitated  AFRL’s  ability  to 
identify  the  weaknesses  of  the  model.  This  work  resulted  in  several  modifications  to  the  IFM 
model. 

2.2  The  Coupled  Ionosphere  Thermosphere  Electrodynamic  Forecast  Model  (CITEFM) 

In  the  current  year,  work  was  initiated  to  incorporate  the  CITEFM  into  applications  at  AFRL. 
The  CITEFM  model  is  a  first-principles  model  that  combines  ionospheric,  thermospheric  and 
dynamo  calculations  to  produce  global  simulations  given  a  set  of  geophysical  conditions.  In  the 
current  year,  we  installed  and  utilized  the  CITEFM  model  to  generate  calculations  of  TEC, 
NmF2,  density  profiles  and  meridional  and  zonal  winds.  These  calculations  were  produced  for  a 


1 


variety  of  geophysical  conditions.  Specialized  software  was  developed  to  graphically  analyze 
and  animate  the  CITEFM  output.  This  work  was  useful  in  detecting  anomalous  behavior  in  the 
model.  As  a  result  of  these  efforts,  the  model  was  revised  to  produce  a  more  accurate 
representation  of  ionospheric  parameters. 

The  work  done  this  year  using  the  IFM  and  CITEFM  was  useful  in  the  initial  plans  of  the 
Communications  /Navigation  Outage  Forecast  System  (C/NOFS)  science  plan.  The  models  were 
used  to  define  and  display  the  current  level  of  knowledge  that  is  built  into  the  specification 
models  used  at  AFRL.  C/NOFS  is  a  science  plan  to  forecast  the  presence  of  equatorial 
ionospheric  irregularities  that  may  adversely  affect  communication  and  navigation  satellite-based 
systems.  It  is  comprised  of  a  satellite  that  will  be  equipped  with  a  variety  of  sensors  to  measure 
electron  densities;  ion  and  electron  temperatures;  magnetic  fields,  neutral  winds,  ionospheric 
scintillations;  and  TEC  along  the  lines  of  sight  between  the  C/NOFS  and  GPS  satellites.  The 
C/NOFS  satellite  will  be  launched  in  early  2004. 

2.3  The  Alfven-Falthammar  Formula  for  the  Parallel  E-Field  and  Its  Analogue  in 
Downward  Auroral-Current  Regions 

In  the  current  year,  efforts  were  expanded  to  explain  the  Birkeland  current  system.  The 
Birkeland  currents  are  a  system  of  upward  and  downward  magnetic  field-aligned  electrical 
currents  in  the  auroral  zone  that  flow  between  the  magnetosphere  and  the  ionosphere  where 
particle  number  flux,  momentum  and  energy  flux  are  exchanged.  The  analysis  used  led  to  a  set 
of  differential  equations.  Several  different  sets  of  differential  equations  were  analyzed.  Each  set 
of  equations  was  analyzed  in  two  steps.  The  first  step  involved  using  symbolic  manipulation  to 
triangularize  the  set  of  equations  to  solve.  The  final  step  of  the  analysis  is  a  Fortran  code,  which 
solves  the  set  of  differential  equations  with  given  boundaries. 

These  formulas  give  new  insight  into  the  physics  of  the  Birkeland  current  system.  It  was  found 
that  for  downward  auroral-current  regions,  the  velocity-space  anisotropy  in  the  ion  distribution 
function  is  such  that  half  the  ion  energy  perpendicular  to  the  magnetic  field  is  greater  than  the 
ion  energy  parallel  to  the  magnetic  field.  This  is  just  opposite  to  the  case  for  upward  auroral- 
current  regions. 

Details  of  this  work  together  with  the  results  were  published  in  the  IEEE  Transactions  on  Plasma 
Science.  A  reprint  of  this  article  is  attached  to  this  annual  report  (see  Appendix). 

2.4  PRISM  Validations 

PRISM  is  a  model  developed  by  AFRL  that  specifies  electron  density  of  the  ionosphere  on  a 
global  and/or  regional  scale.  PRISM  is  capable  of  ingesting  data  from  three  types  of  sensors:  (1) 
ionospheric  sounders;  (2)  TEC  sensors;  (3)  in  situ  sensors  on  board  Defense  Meteorological 
Satellite  Program  (DMSP).  PRISM  then  adjusts  its  climatological  map  to  agree  with  the  data 
taken  at  the  various  geographic  locations. 

In  the  current  year,  we  have  participated  in  continued  efforts  to  determine  the  accuracy  of 
PRISM  calculations.  The  current  efforts  included  defining  the  accuracy  dependence  with  respect 
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to  the  number  of  stations  used  to  drive  the  model.  This  effort  required  collecting  and  analyzing 
data  from  approximately  90  different  geographic  locations.  Our  efforts  were  primarily  devoted 
to  the  analysis  of  the  GPS  TEC  used  in  this  work.  The  data  was  used  to  define  how  well  the 
model  calculations  improve  with  increased  data. 

This  work  was  presented  at  the  AIAA  Space  2000  Conference,  Long  Beach  CA,  September 

2000. 

2.5  The  Great  Magnetic  Storm  of  July  2000 

On  July  14th,  a  coronal  mass  ejection  left  the  Sun  at  a  speed  of  1800  km/s.  Its  effects  on  Earth 
were  noted  on  July  15th  in  the  form  of  an  extreme  geomagnetic  storm  with  Kp  values  of  9 
recorded  from  approximately  15  through  21  UT.  A  storm  of  this  size  can  have  a  significant 
effect  on  TEC  and  consequently  on  sophisticated  space-based  technology  systems,  such  as 
communication  and  navigation.  In  response  to  this  storm,  we  obtained  and  analyzed  worldwide 
GPS  TEC  data  to  determine  the  ionospheric  behavior  during  this  great  storm.  The  most 
significant  findings  were  the  unprecedented  magnitudes  and  gradients  measured  in  TEC  in  the 
southeast  region  of  the  United  States. 


Figure  1.  Equivalent  Vertical  Range  Delay  Measured  from  Miami  During  the  Great 
Magnetic  Storm  of  July  2000. 

This  is  illustrated  in  Figure  1,  where  GPS  measurements  of  equivalent  vertical  range  delay 
recorded  along  multiple  lines  of  sight  from  Miami,  FL  are  plotted  versus  Universal  Time  (UT). 
Note  that  the  range  delays  exceed  40  meters  (equivalent  to  240  TEC  units)  near  23UT.  Also  note 
the  large  spatial  gradients  observed  along  different  lines  of  sight  from  the  Miami  receiver.  All  of 
these  effects  were  unprecedented  features  of  this  mid-latitude  location.  This  work  together 
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additional  storm  analysis  efforts  was  useful  in  defining  mid-latitude  response  to  severe 
geomagnetic  storm  activity. 

2.6  GPS  Position  Error  Maps 

In  the  last  two  years,  we  worked  to  develop  software  to  investigate  navigational  position  errors 
caused  by  inaccuracies  in  the  standard  Ionospheric  Corrections  Algorithm  (ICA).  The  work 
resulted  in  a  model  that  estimates  position  errors  for  the  single  frequency  user  and  displays  these 
potential  errors  on  a  worldwide  graphics  display.  Near  the  end  of  the  last  year,  we  also  worked  to 
incorporate  this  software  into  the  Operational  Space  Environment  Network  Display  (OP-SEND). 
The  OP-SEND  project  is  a  graphic-based  space  environment  forecasting  tool  that  was  made 
operational  at  the  55th  Space  Weather  Squadron  (SWXS).  As  a  result  of  our  efforts  in  the  design 
of  OP-SEND  together  with  our  participation  in  the  transition  of  the  OP-SEND  product  to  SWXS, 
we  were  included  in  the  team  that  was  awarded  the  Annual  Air  Force  Merewether  Award  for  the 
year  2000.  This  award  was  in  the  category  of  Most  Significant  Technical  Contribution. 

OP-SEND  includes  maps  with  near  real-time  estimates  on: 

•  UHF  SATCOM  Scintillation. 

•  HF  Illumination. 

•  Auroral  Clutter  Boundaries. 

•  Scintillation  Predictions. 

•  GPS  Single-Frequency  Position  Errors. 

In  the  current  year,  we  have  also  continued  to  augment  and  support  the  GPS  Single-Frequency 
Position  Errors  component  of  OP-SEND  by  designing  a  comprehensive  validation  effort  to  be 
conducted  over  the  next  two  years  of  this  contract. 

This  work  resulted  in  a  presentation  at  the  American  Geophysical  Union  meeting  held  in  San 
Francisco,  CA  in  December  2000. 

2.7  Ultra-Violet  (UV)  PRISM 

In  the  current  year,  work  was  initiated  to  develop  a  software  product  that  will  be  used  to  integrate 
data  from  various  sensors  into  the  existing  technology  of  the  PRISM  model. 

Two  UV  sensors  are  planned  to  fly  on  the  DMSP  Block  5D-3  satellites.  Both  are  remote-sensing 
instruments  that  measure  ultraviolet  emissions  emitted  by  the  Earth’s  upper  atmosphere.  One 
instrument  is  the  Special  Sensor  Ultraviolet  Limb  Imager  (SSULI),  which  scans  the  limb  of  the 
Earth  in  the  orbital  plane.  It  measures  the  UV  emissions  in  the  interval  from  800  to  1700  A.  The 
other  instrument  is  the  Special  Sensor  Ultraviolet  Spectrographic  Imager  (SSUSI),  that 
measures  emissions  in  five  different  wavelength  bands.  The  resultant  images  cover  the  visible 
Earth  disk  from  horizon  to  horizon  and  the  anti-sunward  limb  up  to  an  altitude  of  approximately 
520  km. 
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These  sensors  have  the  potential  to  significantly  improve  ionospheric  remote  sensing.  To  take 
advantage  of  this  opportunity,  the  data  from  these  new  sensors  must  be  able  to  be  used  by 
operational  products.  Modifying  PRISM  to  accept  new  data  types  is  a  difficult  and  time- 
consuming  task.  Therefore,  we  have  initiated  development  of  software  that  will  accept  the 
DMSP  UV  sensor  data,  extract  the  relevant  data  and  convert  it  to  a  data  type  currently  accepted 
by  PRISM.  This  software  will  be  called  UV  PRISM.  Our  initial  work  indicated  that  the  data 
array  sizes  of  PRISM  were  too  small.  This  resulted  in  a  modified  PRISM  model. 

In  the  current  year,  we  also  simulated  the  UV  measurements  by  using  an  authentic  DMSP  track, 
the  International  Reference  Ionosphere  Model  (IRI)  and  the  Parameterized  Ionosphere  Model 
(PIM).  We  utilized  the  simulated  measurements  as  input  to  PRISM  and  to  develop  graphics  that 
plot  the  results.  These  efforts  allowed  us  to  further  understand  the  effect  that  the  UV 
measurements  will  have  on  PRISM  and  to  determine  how  to  best  utilize  this  new  source  of 
information. 

This  work  is  expected  to  fully  develop  in  the  next  two  years  of  this  contract. 


3.  PRESENTATIONS 

The  work  described  in  this  report  resulted  in  three  oral  presentations  at  various  scientific 
meetings: 

Bishop,  G.,  Quigley,  S.,  Groves,  K.,  Bulled,  T.,  Doherty,  P.,  Cook,  C.,  Citrone,  P.,  Scro,  K.  and 
Wilkes,  R.,  “OpSend:  Mission-Tailored  Graphical  Products  for  DoD  Warfighters,”  presented  at 
the  AGU  Fall  Meeting,  San  Francisco  CA,  December  2000. 

Pulliam,  R.,  Borer,  W.S.,  Decker,  D.T.,  and  Doherty,  P.H.,  "Operational  Ionosphere  Model 
Validations  Study,”  presented  at  the  AIAA  Space  2000  Conference,  Long  Beach  CA,  September 
2000. 

Pulliam,  R.,  Borer,  W.S.,  Decker,  D.T.,  and  Doherty,  P.H.,  "Evaluation  of  Metrics  for  Ambient 
Ionospheric  Specification,”  presented  at  the  AGU  Fall  Meeting,  San  Francisco  CA,  December 
2000. 


4.  JOURNAL  ARTICLES 

Jasperse,  J.R.  and  Grossbard,  N.J.,  "The  Alfven-Falthammer  Formula  for  the  Parallel  E-Field  and 
its  Analogue  in  Downward  Auroral-Current  Regions,”  IEEE  Transactions  on  Plasma  Science, 
Vol.  28,  No.  6,  December  2000. 
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5.  AWARDS 


Patricia  H.  Doherty  and  Susan  H.  Delay  were  part  of  the  team  that  was  awarded  The  Annual  Air 
Force  Merewether  Award  for  the  Year  2000  in  the  category  of  Most  Significant  Technical 
Contribution. 
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The  Alfven-Falthammar  Formula  for  the 
Parallel  .E-Field  and  its  Analogue  in  Downward 
Auroral-Current  Regions 

John  R.  Jasperse  and  Neil  J.  Grossbard,  Member,  IEEE 


Abstract — In  this  paper,  we  give  an  alternative  derivation  of  the 
Alfvln-Falthammar  formula  for  a  positive  parallel  JS-field  in  up¬ 
ward  auroral -current  regions  and  its  analogue  for  a  negative  par¬ 
allel  J£-fieId  in  downward  auroral-current  regions.  These  formulas 
give  new  insight  into  the  physics  of  the  Birkeland  current  system. 
We  find  that  for  downward  auroral-current  regions,  the  velocity- 
space  anisotropy  in  the  ion  distribution  function  is  such  that  half 
the  ion  energy  perpendicular  to  the  magnetic  field  is  greater  than 
the  ion  energy  parallel  to  the  magnetic  field.  This  is  just  opposite 
to  the  case  for  upward  auroral-current  regions.  These  results  are 
compared  to  recent  particle-in-cell  simulations  and  FREJA  satel¬ 
lite  data. 

Index  Terms — Alfven-Falthammar  formula,  Birkeland  current 
system,  parallel  ^-fields  in  space  plasmas. 


I.  Introduction 

THIS  paper  deals  with  the  theory  of  the  Birkeland  current 
system  which  was  discovered  nearly  100  years  ago  and  re¬ 
mains  an  active  area  of  space  physics  research  to  this  day.  The 
Birkeland  currents  are  a  system  of  upward  and  downward  mag¬ 
netic  field-aligned  electrical  currents  in  the  auroral  zone  that 
flow  between  the  magnetosphere  and  the  ionosphere  where  par¬ 
ticle  number  flux,  momentum,  and  energy  flux  are  exchanged. 
In  our  view,  there  is  no  problem  more  important  to  the  under¬ 
standing  of  space  weather  at  high  latitudes  than  to  understand 
how  the  Birkeland  current  system  works.  For  a  study  of  the  sta¬ 
tistical  properties  of  the  Birkeland  currents  using  TRIAD  satel¬ 
lite  data,  see  [1,  Fig.  13].  For  a  recent  case  study  using  FAST 
satellite  data,  see  [2,  Fig.  3]. 

One  of  the  purposes  of  this  paper  is  to  give  an  alternative 
derivation  of  the  Alfven-Falthammar  ( A-F)  formula  for  the  par¬ 
allel  £-field  (E||)  in  upward  auroral-current  regions.  There  are 
several  ways  to  derive  this  result  [3]— [6],  The  formulas  for  E\\ 
are  the  same  but  the  derivations  are  not  equivalent  since  different 
assumptions,  closure  approximations,  and  boundary  conditions 
are  used.  Our  derivation  begins  with  the  Vlasov  equations  and, 
as  we  will  show,  not  only  gives  new  insight  into  the  physics  of 
that  problem,  but  may  be  generalized  in  a  straightforward  way 
to  study  downward  auroral-current  regions. 
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energetic 

precipitating 


potential  decrease  from  s,  to  s2  ($2  -  <  0) 

downflowing  ’’inverted-V*  and  field-aligned  electrons 
upfbwing  ion  beams 
transverse  ion  heating 

electric-and  magnetic-field  turbulence  and  ion  cydotron  waves 
Fig.  1 .  Some  properties  of  upward  auroral -current  regions. 

The  primary  purpose  of  this  paper,  therefore,  is  to  give  the 
analogue  of  the  A-F  formula  for  E\\  for  downward  auroral-cur- 
rent  regions.  The  derivation  begins  with  a  velocity-space  diffu¬ 
sion  model  for  the  problem,  and  the  formula  that  we  give  for  E\\ 
also  gives  new  insight  into  the  physics  of  that  problem. 

There  is  a  vast  literature  on  upward  auroral-current  regions. 
For  a  review  of  rocket  and  satellite  measurements  before  1980, 
see  [7],  For  a  sample  of  more  recent  satellite  measurements, 
see  [8]  on  FREJA  results  and  FAST  results  in  {9].  For  example, 
the  properties  of  upward  auroral-current  regions  have  been  ob¬ 
served  by  the  ISIS-1  [10],  S3-3  [7],  NOAA-6  [1 1],  DE-1  [12], 
DE-2  [13],  FREJA  [14],  POLAR  [15],  and  FAST  [16]  satel¬ 
lites,  to  name  a  few.  From  these  studies,  we  see  that  upward 
auroral-current  regions  are  often  characterized  by 

1)  converging  electrostatic  shocks,  upward  pointing  E\\\ 

2)  downflowing  “inverted-V”  and  field-aligned  electrons; 

3)  large-scale  density  cavities; 

4)  upflowing  ion  beams  and  transverse  ion  heating; 

5)  electric-  and  magnetic- field  turbulence  and  ion  cyclotron 
waves; 

6)  large  fluxes  of  earthward-streaming  energetic  (>  few 
kilovolts)  magnetospheric  particles. 

Some  of  these  properties  are  illustrated  in  Fig.  1  (not  to  scale). 
The  theoretical  literature  before  1993  is  reviewed  in  [17).  For  a 
more  recent  review,  see  (18,  Introduction  and  Bibliography]. 

The  properties  of  downward  auroral-current  regions  are  less- 
explored.  They  have  been  observed  by  the  ISIS-2  [19],  DE-1 
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t  H7^ 


E.<0 


.  small  fluxes  of 
/  magnetospheric  parides 


upfiowing  field-aligned  electrons 
upfiowing  energetic  ion  conics 


potential  increase  from  s,  to  s2  (+2  -  >  0) 

upfiowing  field-aligned  electrons 
intense  transverse  ion  beating 
intense  ELF  and  VLF  electric-field  turbulence 
fast  bipolar  solitary  structures 

Fig.  2.  Some  properties  of  downward  auroral -current  regions. 

[20],  S3-3  [21],  FREJA  [22],  FAST  [23],  [24],  and  other  satel¬ 
lites.  These  downward  currents  often  occur  in  the  auroral  zone 
adjacent  to  “inverted- V”  structures  [2],  the  dayside  cusp-cleft 
region  [21],  and  in  association  with  the  “black  aurora”  [25]. 
From  these  studies,  we  see  that  downward  auroral -cunent  re¬ 
gions  are  often  characterized  by 

1)  diverging  electrostatic  shocks,  downward  pointing  Ey; 

2)  upfiowing  low  energy  (<  few  kilovolts)  Field-aligned 
electrons; 

3)  intense  transverse  ion  heating; 

4)  small-scale  density  cavities; 

5)  intense  ELF  and  VLF  electric-field  turbulence  and  fast 
bipolar  solitary  structures; 

6)  small  fluxes  of  energetic  (>  few  kilovolts)  magneto- 
spheric  particles. 

Some  of  these  properties  are  illustrated  in  Fig.  2  (not  to  scale). 
Some  theoretical  work  has  been  carried  out  for  downward  au¬ 
roral-current  regions  using  the  multimoment  fluid  equations  for 
the  current-driven  ion  cyclotron  instability  [26];  particle-in-cell 
(PIC)  simulations  in  the  electrostatic  approximation  [27],  [28]; 
and  a  velocity -space  diffusion  model  for  the  problem  [29],  [30]. 

In  Section  II,  we  give  the  alternative  derivation  of  the  A-F 
formula  for  a  positive  Ey  in  an  upward  auroral-current  region, 
and,  in  Section  III,  its  analogue  for  a  negative  Ey  in  a  downward 
auroral-current  region.  In  Section  IV,  we  discuss  these  results. 


presence  of  low-energy  (or  ambient)  ionospheric  particles, 
and  wave-particle  interactions.  The  simplest  way  to  model 
an  auroral  arc  that  wc  can  think  of  is  to  neglect  the  ambient 
ionospheric  particles  and  wave-panicle  interactions,  and  to 
consider  only  the  injected  magnetospheric  particles.  In  his 
review,  Borovsky  [17]  discusses  12  mechanisms  by  which 
the  auroral  plasma  can  sustain  an  Ey,  thereby  accelerating 
aurora!  particles,  and  also  ten  drivers  (or  generators)  of  the 
mechanisms.  As  we  will  see  in  this  section,  the  A-F  formula 
for  a  positive  Ey  in  upward  auroral-current  regions  results 
from  the  following  simple  model  for  a  quasi-stable  auroral  arc 
I)  the  driver  is  the  earthward-streaming  energetic  magneto¬ 
spheric  particles  and  2)  the  mechanism  that  sustains  Ey  is  the 
differential  velocity-space  anisotropy  between  the  ions  and  the 
electrons  in  the  dipolar  magnetic  field. 

In  this  paper,  we  assume  that  the  gyrotropic  approximation  is 
valid  and  seek  a  steady-state  solution.  The  Vlasov  equations  in 
one-spatial  and  two-velocity  dimensions  are 


f  a 

\v"Ts  + 


d 

dv\\ 


_L  — 

2  B  ds 


d 


and  the  quasi-neutrality  (Q-N)  condition  is 

nc  =  ^  np 
3 


(1) 

(2) 


where  fi  is  summed  over  the  ions.  Here,  py  and  Ey  are  the  par¬ 
allel  components  of  the  gravitational  acceleration  and  the  elec¬ 
trostatic  field,  respectively;  B  is  the  geomagnetic  field;  /Q  = 
fQ(s, t;x,  vjj)  is  the  gyrotropically-averaged  distribution  func¬ 
tion  for  each  species  a;  and  5  is  the  distance  from  the  earth  along 
the  geomagnetic  flux  tube.  Also,  fQ  is  normalized  to  the  particle 
density,  ntt,  and  the  other  quantities  have  their  usual  meanings. 

In  (2),  we  have  replaced  Poisson’s  equation  with  the  Q-N  con¬ 
dition.  This  greatly  simplifies  the  problem.  The  standard  argu¬ 
ment  for  making  this  replacement  is  that  for  length  scales  which 
are  long  compared  to  the  Debye  length,  the  two  are  equiva¬ 
lent.  However,  this  is  not  at  all  obvious  since,  in  making  this 
replacement,  we  have  completely  changed  the  nonlinear  char¬ 
acter  and  boundary  conditions  of  the  problem.  Much  has  been 
written  about  this  issue  [4]-[6],  [31  ]— [33]. 


II.  Parallel  E-Field  Due  to  Injected  Magnetospheric 
Particles  in  Upward  Auroral-Current  Regions 

The  purpose  of  this  Section  is  to  give  an  alternative  deriva¬ 
tion  of  the  AIfv£n-Falthammar  (A-F)  formula  for  the  parallel 
E-field  (Ey)  in  upward  auroral-current  regions.  The  original 
derivation  is  given  in  [3].  Although  we  obtain  the  same  formula 
for  Ey,  the  method  we  use  is  different.  Our  reason  for  presenting 
an  alternative  derivation  here  is  that  the  method  we  use  may  be 
easily  generalized  and  applied  to  downward  auroral -current  re¬ 
gions,  as  we  will  see  in  Section  III. 

A.  Vlasov  Equations  and  the  Quasi-Neutrality  Condition 

Upward  auroral -current  regions  are  often  characterized  by 
earth  ward- streaming  energetic  magnetospheric  particles,  the 


B.  A-F  Solution  for  Ey  and  Related  Quantities  from  the 
Multimoment  Fluid  Equations 

We  now  proceed  to  derive  the  multimoment  fluid  equations 
by  multiplying  (1)  by  v"  vjj  and  integrating  over  all  velocity 
space.  We  use  the  following  notation: 


/oo  roc 

dv ||  /  i>x,V||) 

•  oo  Jo 

A  =  -B~l  dB/ds  =  A~l  dA/ds 


(3) 

(4) 


where  A  is  the  cross  sectional  area  of  the  flux  tube.  Here  n  and 
i  are  positive  integers  or  zero.  We  obtain  for  all  n  and  l  =  0 

d/ds(n .  1)„  +  A(  1  +  n/2)(n,  l)n  =  0  (5) 
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and  for  all  n  and  f  >  I 

d/dfi(n.f  4-  I)tt  +  A(i  +  n/2){n,f  +  l)„ 

+  f{{fl,j™0)(d/ds$',u)(n. .(  -  1)„ 

-  (A/2){n  +  2.t  -  l)n  }  =  0.  (6) 

Here,  $0|(  =  <j>\\  +  (ma/qa)<pc g,  where  <t> t  is  the  electrostatic 
potential  and  <t>G\\  is  the  gravitational  potential.  The  multimo- 
ment  fluid  equations  for  n  =  0,  f.  =  0;  n  =  0,  l  =  1;  n  =  0, 


(  =  2;  and  «  =  2,  l  =  0  are 

d/ds  nauQ  +  Anaua  =  0  (7) 

rf/d-s  n„v)Qii  +  An„(wa||  -  waL/2) 

+  (Va/2 )na  d/ds$a\\  =  0  (8) 

d/ds  7t„Q„||  +  J4n„(Q0||  -  Qa±) 

+  qan0u0  d/ds$a ||  =  0  (9) 

d/ds  naQa j.  +  2i47i„Qox  =  0-  00) 


Here,  we  have  introduced  the  standard  definitions 

(0, 0)o  =  71q  i  (0.  l)a  “  naua|  mQ(2,0)a/2  = 

ma(0,2)a/2  =  n„wQ||;  ma(2,l>a/2  =  naXQax;  and 
ma(0,3)a/2  =  naQajj.  The  Q$  denote  the  energy  flux  and 
should  not  be  confused  with  the  heat  flux.  Also,  by  definition, 
wol  =  T„_l  and 

walj  =Toll/2  +  maiil/2 

where  the  temperatures  are  expressed  in  energy  units.  Adding 
(9)  and  ( 10)  together  and  using  (7),  we  may  obtain  an  alternative 
equation  for  (9).  Equations  (7H10)  may  then  be  rewritten  as 

B  d/ds(naua/B)  =  0  (11) 

d/ds(nawa\\)  -  B”1  dB/ds  na(ti/aj|  -  wa±/2) 

4*  nQ(qa/ 2)  d/ds  $a|j  =  0  (12) 

B  d/ds{[naQQ ||  +  HaQa±  4“  ^a^o^o^aJl]/^}  ^  0  (13) 

B2  d/ds(nQQa±/B2)  =  0.  (14) 

Equations  (1 1H14)  are  the  standard  conservation  laws  for  the 
number  flux,  momentum,  total  energy  flux,  and  perpendicular 
energy  flux,  respectively. 

Formal  expressions  for  E\\  may  be  found  from  (12)  and  (13) 
by  simply  rearranging  the  terms.  Neglecting  gravity,  we  see 
from  momentum  conservation  (12)  that 

£(|  =(2 /qa){-B-'  dB/ds{w^  -W.J./2) 

+  n2l  d/ds  naiya|[}  (15) 

and  from  total  energy  flux  conservation  (13)  with  the  help  of 
(11)  that 

E\\  =  (B/qQnaua)  d/ds[(naQ q)!  4-  naQa±)/B].  (16) 

Equations  (15)  and  (16)  are  2a  expressions  for  E y  in  terms  of 
the  moments  of  fQ.  We  call  them  self-consistent  formulas  for 
E|l  since,  once  the  explicit  solution  for  the  moments  of  fQ  are 
found,  then  (15)  and  (16)  must  be  satisfied. 

Equations  (1 1M14)  and  (2)  are  4a  +  1  equations  for  6a  + 
1  unknowns  including  the  electrostatic  potential  </>\ j.  In  order 
to  solve  them,  we  must  introduce  closure  approximations.  The 
standard  procedure  for  closing  the  multimoment  fluid  equations 
is  to  use  the  gyrotropic  limit  of  the  Maxwellian-based  5-,  8-,  10-, 
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or  13-moment  sets,  or  the  bi-Maxwellian-based  6-  or  16-mo¬ 
ment  sets.  For  a  discussion  of  this  procedure,  see  |34|.  Specific 
solutions  have  also  been  worked  out  for  a  variety  ol  plasma- How 
conditions  |35|.  [36|.  Instead  ol  using  the  standard  procedure, 
we  use  a  simpler  closure  approximation  which,  as  we  shall  see 
later  in  this  section,  places  a  restriction  on  the  type  of  flow  con¬ 


ditions  to  which  it  may  be  applied.  The  closure  wc  use  is 

Qa±=U  aV>a±  (17) 

Qa||— U®  Wo||*  08) 

Using  (I  I),  (17),  and  (18),  we  find  that  the  multimoment  fluid 
equations  become 

d/ds(nQuQ  /  B)  =  0  (19) 

d/ds  n0iuQ(j  -  B~}  dB/ds  n0(wQ||  -  w«J2) 

+  »UW2)*t/rf«*«„=0  (20) 

d/ds(wQ ||  +  wa±  4*  qQ^a\\)  =  0  (21) 

d/ds(wa±/B)  =  0  (22) 

ne  ng,  0  =  ions.  (23) 

e 


A  few  remarks  about  (I9H23)  are  in  order  here.  First,  they  are 
a  set  of  4a  4“  1  nonlinear,  first-order  ordinary  differentia!  equa¬ 
tions  for  4a  +  1  unknowns  subject  to  boundary  conditions  im¬ 
posed  at  either  end  of  the  flux  tube.  Second,  the  closure  approx¬ 
imation  that  we  have  used  has,  in  effect,  replaced  the  conser¬ 
vation  of  total  energy  flux  (13)  and  the  conservation  of  perpen¬ 
dicular  energy  flux  (14),  both  of  which  are  rigorous,  by  conser¬ 
vation  relations  for  the  total  energy  (21)  and  the  perpendicular 
energy  (22).  These  are  just  the  conservation  relations  that  we 
obtain  from  single-particle  orbit  theory.  However,  the  concept 
of  parallel  and  perpendicular  temperature  has  been  preserved  in 
(I9M23)  as  wa±  =  Tq_l  and  u;a||  =  Ta||/24*mau2/2.  Third, 
when  anomalous  transport  effects  are  neglected,  (19)  and  (20) 
are  identical  to  the  steady-state  versions  of  (1)  and  (2)  in  [37], 
where  the  authors  have  used  the  gyrotropic  limit  of  the  16-mo¬ 
ment  set  to  study  the  polar  wind.  However,  (21)  and  (22)  differ 
from  the  higher  order  moment  equations  in  [37],  as  the  closure 
approximation  is  different. 

Now,  we  proceed  to  solve  (19M23)  in  closed  form  subject  to 
appropriate  boundary  conditions.  Consider  s\  <  s  <  where 
Si  is  in  the  ionosphere  and  $2  is  in  the  magnetosphere.  Equation 
(20)  may  be  further  simplified.  We  may  eliminate  $a\\  from  (8) 
and  (9)  and  use  closure  to  obtain 

d/ds  nQuftw0||  -  2u0  d/ds  nawa\\  -  AnQuQwQ\\  =  0. 

Using  (7),  we  may  further  simplify  to  obtain 

2iuq||  dna/ds  4-  2Anawa^  4*  na  dwQ\\/ds  =  0. 

The  left-hand  side  of  this  expression  is  proportional 
to  d/ds(n2wa\\/B2)  which  is  in  turn  proportional  lo 
d/ds(wQft/ul)t  so  (20)  may  be  replaced  by 

d/ds(wal]/u2Q )  =  0.  (24) 

The  hidden  meaning  of  the  closure  approximation  given  by  ( 1 8) 
is  now  revealed.  If  the  ratio  (TQ||/2  4*  rnnu2a/2)l{rnau2Ct/2) 
does  not  change  by  very  much  on  the  flux  tube,  then  (18)  is 
a  reasonable  approximation.  This  requires  that  Tft||  and  m.0u2a 
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change  proportionally  from  *j  to  *2  so  l^a!  a^ovc  rat'°  *s 
nearly  constant.  If  we  examine  the  experimental  data  and  find 
that  this  approximation  is  valid,  then  ( 18)  is  a  good  closure  ap¬ 
proximation.  However,  we  note  that  this  closure  may  not  apply 
to  problems  where  the  plasma  flow  evolves  from  subsonic  to 
supersonic.  We  now  integrate  (22),  (21 ),  (24),  and  (19)  from  .s 
to  ,s2,  in  that  order,  to  find  wol%  ti/a|j,  u0  and  n0  in  terms  of 
4> ||  and  the  boundary  values  at  s2.  We  then  specialize  to  an  elec¬ 
tron-ion  plasma  and  impose  the  Q-N  condition  to  solve  for  <j> ||. 
Neglecting  gravity,  the  solution  is 


Wal(»)  =  UlQj_2[B(s)/B2] 

1«a||(*)  =™o||2  +  ti70 j_2(l  ~  B(s)/B2] 

+  </a[$||2  -  0||M] 

ua(s)  =Ua2[li;Q||(.s)/l/7a||2]1/2 

na(s)  =Uft2(5(5)/^2)[^ol|2/w;a||(«)]1/2 


>e||2^tX2l  \B(s)  _  x 

,Je\\2)  J  #2 


(25) 

(26) 

(27) 

(28) 

(29) 


e{wm  +  we 

where  the  subscript  2  denotes  the  value  of  the  quantity  at  s2  and 
e  is  the  magnitude  of  the  electronic  charge.  From  (29),  we  see 
that  £||  is 


En  =  -  d^/ds  =  -K  dB/ds  (30) 


r 1 ) 

U’»||2UV-L2  ~  t0e||2tl/,X2' 

\eB2J 

W,||2  +  We  ||2 

Equation  (30)  is  the  A-F  formula  for  E\\  for  upward  auroral-cur¬ 
rent  regions  and  is  identical  to  ( 1 2)  of  [3,  Sec.  5.1 .3],  except  that 
™Qj|  has  both  a  thermal  part  and  a  drift  part,  as  wa\\  -  Ta  j|/2  + 
mau2a  /2.  It  was  derived  in  [3]  by  a  different  method,  i.e.,  single- 
particle  orbit  theory,  and  applied  for  different  boundary  condi¬ 
tions.  The  boundary  conditions  used  in  [3]  are  that  ionospheric 
ions  are  emitted  at  the  lower  boundary  and  magnetosperic  elec¬ 
trons  are  injected  at  the  upper  boundary.  For  that  problem,  tin 
and  hence  U&  are  positive,  and  ue2  is  negative.  In  this  paper, 
we  have  applied  the  A-F  formula  to  the  problem  where  the 
ions  and  electrons  are  of  magnetospheric  origin  and  both  are 
injected  at  the  upper  boundary,  hence  Ui2  and  uc2  are  negative. 
Using  (25),  (26),  and  (29),  we  may  also  show  that  K  is  an  in¬ 
variant.  Equation  (30)  is  also  equivalent  to  (28)  obtained  in  [6] 
by  using  Persson’s  method  of  solution  [4],  [5].  As  pointed  out  in 
[3],  the  differential  particle  anisotropy  is  such  that  [ii7;||2twex2  - 
t/7e|i2u7jj_2]  is  positive  and,  since  dB/ds  is  negative,  E\\  is  up¬ 
ward  pointing. 

A  graph  of  the  solution  given  by  (25M30)  is  shown  in  Fig.  3, 
where  s  measures  the  distance  along  B  from  the  surface  of  the 
earth.  Remember  that  tu0x  =  Tax-  We  have  also  plotted  the 
current  density,  j||,  which  is  defined  as  =  en(ui  —  ue). 
Here  we  have  used  the  superscript  M  to  denote  the  fact  that 
the  particles  are  of  magnetospheric  origin.  In  plotting  the  so¬ 
lution,  we  have  imposed  the  following  conditions.  The  phys¬ 
ical  conditions  are  $i  —  «o  =  1000  km,  s2  —  ao  —  6500 
km,  rrii/rnc  =  250,  B(s)  =  Bo{so/s)3,  B0  =  0.57  Gauss, 
where  s0,  the  radius  of  the  earth,  is  approximated  by  6400  km. 
The  boundary  conditions  at  a2  are  rtj2  =  nr2  =  n2  =  0.50 
cm"3,  ui2  =  —5.6  x  103  km/s,  uc2  =  -6.8  x  103  km/s, 
Ti||2  =  Te|{2  =  TiX  2  =  2  =  1.0  keV,  and  $„2  =  -4.1  x  103 


V  These  parameters  give /||2  =  0.10  /tA/nr,  u>,||2  =  22 A  keV, 
and  «Vj|2  =  0.62  keV  Since  $ j|  is  arbitrary  to  a  constant,  we 
have  adjusted  $|i2  so  that  =  0  for  convenience  in  plotting 
The  mass  and  temperature  ratios  are  unrealistic  and  the  current 
density  is  too  small.  However,  w^2  =  22.4  keV  is  realistic  for 
strong  injection  events  [38].  We  use  these  parameters  in  order  to 
compare  the  A-F  solution  to  Schriver’s  simulation  [  18)  in  Sec¬ 
tion  II-D.  In  this  example,  since  u;ex2  is  comparable  to  ?i^x2 
and  u/jj|2  »  uv}|2»  a  good  approximation  for  $j|  is 

<l>\\{s)  =  (t>\\ 2  +  e_liuex2[7(*)  “  [\  (32) 

where  7(a)  =  B(s)/B2.  For  these  conditions,  the  potential  drop 
from  6*i  to  a2  (<^||2  -  $||i)  scales  as  ujcx2- 

C.  Self -Consistent  Formulas  for  E^ 

Equations  ( 1 5)  and  ( 1 6)  give  four  self-consistent  formulas  for 
£|l  in  an  electron-ion  plasma.  In  order  to  verify  the  validity  of 
these  expressions,  we  use  closure,  substitute  (25>— (29)  into  the 
right-hand  sides  of  (15)  and  (16),  and  find  that  E\\  is  indeed 
given  by  (30). 

Consider  (15)  for  the  ions.  It  is 

E\\  =(2/e){-£"1  dB/ds[w^  -wixl 2] 

+  nf1  d/ds  riiu;l||}.  (33) 

We  wish  to  demonstrate  that  E\\  is  positive  by  using  (33)  and  the 
results  given  in  Fig.  3  instead  of  (30).  To  do  this,  we  note  that 
since  —  B~l  dB/ds  is  positive  and  luqj  >  Wi±/2,  then  the  first 
term  on  the  right-hand  side  of  (33)  is  positive.  The  second  term 
.is  negative,  but  {nf1  d/ds‘  njtUj|||  is  smaller  than  the  first  term 
inside  the  brackets,  so  E\\  is  positive.  This  is  shown  in  Fig.  4, 
where  E^  and  denote  the  first  and  second  terms  on  the 

right-hand  side  of  (33),  respectively,  and  Ejjtota^  is  the  total  E\y 
We  see  that  once  the  solution  is  known,  then  the  fact  that  E\\  is 
positive  can  be  seen  from  the  ion  moments  alone. 

D.  Comparison  to  the  Schriver  Particle-in-Cell  Simulation 

In  this  section,  we  compare  the  A-F  solution  to  the  par- 
ticle-in-cell  (PIC)  simulation  in  the  electrostatic  approximation 
given  by  Schriver  [18].  The  physical  and  boundary  conditions 
that  we  chose  in  Section  II-B  are  essentially  identical  to  run 
number  3  in  [18,  Table  2],  where  Ui/vte  at  s2  is  0.3.  In  both 
cases,  Ui 2  =  -5.6  x  103  km/s  and  w^2  =  22.4  keV.  The 
fact  that  j||2  was  0.1  pAJm2  in  Section  II-B  and  zero  in  the 
simulation  makes  little  difference  in  the  results. 

As  we  saw  in  Section  II-B,  the  A-F  formula  for  E\\  was  ob¬ 
tained  when  the  ionospheric  particles  and  wave-particle  inter¬ 
actions  were  neglected,  the  Q-N  condition  was  used  in  place 
of  Poisson’s  equation,  and  only  the  injected  energetic  magneto¬ 
spheric  particles  were  included.  From  [18,  Fig.  5],  we  see  that 
the  potential  drop  is  —2.7  kV  compared  to  —4.1  kV  from  the 
A-F  formula  shown  in  Fig.  3,  and  that  the  shape  of  $||(s)  is 
concave  downward  (d2<t>\\/ds2  <  0)  instead  of  convex  upward 
(fp^i/ds2  >  0)  in  the  A-F  formula.  To  find  the  polarization 
charge  density,  we  may  substitute  $u  into  Poisson’s  equation. 
We  find  that  the  shape  of  $j|  in  the  simulation  is  consistent  with 
the  correct  polarization  charge  density,  whereas  the  A-F  result  is 
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Fig.  3.  The  A-F  solution  for  injected  magnetospheric  particles  in  an  upward  auroral -current  region  from  (25H30).  The  symbols  are  defined  in  the  text,  and  the 
superscript  A/  denotes  particles  of  magnetospheric  origin. 


not.  Apparently,  including  ionospheric  particles  and  wave-par¬ 
ticle  interactions,  and  using  Poisson’s  equation  instead  of  the 
Q-N  condition,  have  modified  the  shape  and  magnitude  of  the 
potential  drop. 

However,  it  appears  that  the  A-F  formula  gives  the  essen¬ 
tial  physics  of  the  problem  correctly,  i.e.,  that  for  earthward- 
streaming  energetic  magnetospheric  particles,  it  is  the  differen¬ 
tial  velocity-space  anisotropy  between  the  ions  and  the  electrons 
that  produces  the  potential  drop  and  the  upward  pointing  E\\. 

III.  Parallel  E-Field  Consistent  with  the 
Turbulent  Heating  of  Ionospheric  Ions  in  Downward 
Auroral-Current  Regions 

The  purpose  of  this  section  is  to  give  the  analogue  of  the  A-F 
formula  for  Ejj  in  downward  auroral -current  regions.  For  up¬ 


ward  currents,  E\\  is  positive,  as  we  saw  in  Section  II;  but  for 
downward  currents,  Ej|  is  negative,  as  we  will  see  in  this  sec¬ 
tion.  We  obtain  this  result  by  generalizing  the  method  presented 
in  Section  II. 

A.  Kinetic  Diffusion  Equations  and  the  Quasi- Neutrality 
Condition 

Downward  auroral-current  regions  are  often  characterized  by 
energetic  ion  conics  and  upflowing,  field-aligned  electrons  of 
ionospheric  origin;  intense  broadband  ELF  (BBELF)  electric- 
field  turbulence;  and  small  fluxes  of  earthward-streaming  mag¬ 
netospheric  particles.  The  simplest  model  for  downward  au¬ 
roral-current  regions  that  we  can  think  of  is  neglecting  the  mag¬ 
netospheric  particles  and  considering  only  the  ionospheric  par¬ 
ticles  and  the  BBELF  turbulence.  In  this  paper,  we  also  assume 
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Fig.  4.  The  solution  for  E\\  in  an  upward  auroral -current  region  from  the  self-consistent  A-F  formula  for  E\\  given  by  (33).  Here,  Ej,1 1  is  the  positive  contribution, 
E[|2)  is  the  negative  contribution,  and  Ej is  the  total. 


that  the  spectrum  of  the  BBELF  turbulence  is  given  experimen¬ 
tally.  How  the  turbulence  is  driven  (or  generated)  is  an  active 
area  of  research  and  beyond  the  scope  of  this  paper.  If  an  ad¬ 
ditional  source  of  ion  heating  exists  on  the  flux  tube  due,  for 
example,  to  the  presence  of  fast  bipolar  solitary  structures  as 
suggested  in  [23]  and  [24],  then  see  Section  IV  for  a  brief  dis¬ 
cussion  of  how  this  effect  may  be  incorporated  into  the  method 
of  solution  presented  in  this  paper. 

In  Section  II-C,  we  showed  that  a  positive  Ey  occurs  when 
Wi\\  >  Wi±/2  for  the  ions,  and  the  magnitude  of  the  second 
term  in  (33)  is  smaller  than  the  first  term.  In  Section  m-C,  we 
will  show  that  a  negative  E\\  will  occur  when  the  reverse  ion 
anisotropy  is  true,  i.e.,  when  Wi±/2  >  and  certain  other 
conditions  are  satisfied.  Thus,  in  this  section,  we  propose  the 
following  model  for  a  negative  E\\  in  a  quasi-stable  downward 
auroral-current  region:  1 )  the  driver  of  the  BBELF  turbulence  is 
unspecified,  but  the  spectrum  of  the  turbulence  is  given  experi¬ 
mentally  and  2)  the  self-consistent  E\\  is  sustained  by  a  two-fold 
mechanism  where  the  BBELF  turbulence  resonantly  heats  the 
ions  perpendicular  to  B  which,  in  a  dipolar  magnetic  field,  pro¬ 
duces  a  velocity-space  anisotropy  in  the  ion  distribution  func¬ 
tion  where  wijL! 2  >  w{\\.  Of  course,  the  negative  E\ \  must 
also  be  self-consistent  with  the  upward  acceleration  of  the  elec¬ 
trons,  and  this  relationship  is  also  discussed  in  Sections  III-B 
and  III-D. 

The  steady-state  kinetic  diffusion  equations  in  one-spatial  and 
two- velocity  dimensions  are 


{ 


d_ 

dv\\ 


_L  ^ 

2  B  ds 


(34) 


(£)„ 

and  the  quasi -neutrality  (Q-N)  condition  is 

«*  =  51  n0 

where  / 3  is  summed  over  the  ions.  Equations  (34>-(36)  are  iden¬ 
tical  to  (1)  and  (2)  except  for  the  velocity-space  diffusion  term. 
Here,  Da±  is  the  perpendicular  component  of  the  gyrotropically 
averaged  diffusion  tensor,  which  is  discussed  below.  Equation 
(34)  has  been  used  before  to  develop  a  kinetic  theory  of  the 
“pressure  cooker”  effect  in  downward  auroral-current  regions 
[29],  [30].  For  the  original  reference  on  the  “pressure  cooker” 
effect,  see  [21]. 

As  mentioned  above,  we  assume  that  BBELF  turbulence  ex¬ 
ists  on  the  flux  tube,  is  given  experimentally,  and  that  the  waves 
heat  the  ions  by  cyclotron  resonance  near  the  ion  gyrofrequency. 
Also,  we  do  not  specify  the  driver  of  the  turbulence  nor  do  we  at¬ 
tempt  a  self-consistent  treatment  of  the  particles  and  the  waves. 

In  general,  the  diffusion  tensor  for  the  problem  contains  both 
a  perpendicular  and  a  parallel  part  and  is  velocity  dependent. 
It  may  be  derived  for  electromagnetic  modes  from  quasi-linear 
theory  [39]  or  from  Fokker-Planck  theory  by  assuming  that  the 
particles  undergo  a  random  walk  in  velocity  space  [40].  If  the 
dominant  modes  are  electromagnetic  and  the  ion  energies  are 
less  than  about  1  KeV,  then  k\\V\\  <  ft,  and  D±  »  Dp  where 
is  the  parallel  wave  vector,  ft*  is  the  ion  gyrofrequency,  and 
D|j  is  the  parallel  component  of  the  diffusion  tensor.  For  these 
conditions,  the  turbulent  (anomalous),  resonant  perpendicular 
ion  heating  rate  is 

wo.  =  2  miDix  =  {q} /2mi)S  l(lj  —  Qi)  (37) 


(35) 

(36) 
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where  .SY(U,)  is  the  fraction  ol  the  total  electric-field  spectral 
density  in  the  lelt-hand  circular  polarized  (LHCP)  electromag¬ 
netic  ion  cyclotron  component  and  is  velocity  independent.  This 
result  may  also  he  derived  from  simple  arguments  [4I|.  How¬ 
ever,  if  the  dominant  modes  are  electrostatic,  then  the  diffusion 
tensor  for  the  problem  contains  velocity  dependence,  and  (35) 
and  (37)  arc  approximate. 

For  the  electrons,  Dt ±  =  0,  since  there  is  no  resonant  heating 
at  these  low  frequencies.  However,  preliminary  studies  show 
that  under  certain  conditions,  nonresonant  wave-electron  inter¬ 
actions  are  important.  In  this  paper,  we  neglect  all  turbulent  elec¬ 
tron  heating.  Therefore,  the  approximation  that  we  use  is  that 
the  kinetic  diffusion  equation  is  used  for  the  ion  dynamics  where 
only  the  perpendicular  diffusion  term  is  retained,  and  the  Vlasov 
equation  is  used  for  the  electron  dynamics. 

In  (36),  Poisson's  equation  has  been  replaced  by  the  Q-N  con¬ 
dition.  The  same  remarks  made  in  Section  II-A  about  the  va¬ 
lidity  of  this  replacement  apply  here. 


B.  Numerical  Solution  for  E\\  and  Related  Quantities  from  the 
Multimoment  Fluid  Equations 

We  now  proceed  to  derive  the  multimoment  fluid  equations 
by  multiplying  (34)  by  vjvjj  and  integrating  over  all  velocity 
space.  We  use  the  same  notation  as  in  Section  II-B.  In  addition, 
we  need  to  evaluate  the  following  integrals: 

2*  J ^  <b|«f  jf  dvLvl-^-vLDaX-^-!a. 

Since  DQ±.  depends  only  on  s  in  our  model,  the  integrals  are 
zero  for  n  =  0,  are 

/oc  roo 

/  dv±fa  =  Dax{-l,t)a 

-oc  JO 


for  n  =  1  and  n2Daj_(n  -  2,  i)Q  for  n  >  2.  Using  these  defini¬ 
tions  for  the  bracket  symbols,  the  multimoment  fluid  equations 
for  all  n  and  i  are 


dj  d$(n,  t  +  l)o  +  >4(1  +  n/2)(n,  i  +  1)« 

+  t{{qa/mQ)(d/ds$a\\)(n,  i  -  l)a 
-(A/2)(n  +  2,*-l)a} 

=  n2  Da±{n  —  2,£)0.  (38) 


Using  the  same  definitions  and  procedures  as  in  Section  II,  we 
find  that  the  multimoment  fluid  equations  for  n  =  0,  i  =  0; 
n  =  0,  £  =  1;  n  =  0,  t  =  2;  and  n  =  2,  t  =  0  are 


B  d/ds(nQua/B)  =  0  (39) 

d/ds(nQwa\\)  -  B*1  dB/ds  ntt(iuQ||  -  waX/2) 

+  nQ(qa/2)  d/ds  $q(|  =  0  (40) 

B  d/ds{\naQQ ()  +  nQQa±  +  qQnaua$o]]\/B} 

=  wQ  xnQ  (41) 

B2  d/ds(naQa±/B2)  =  wQ±na.  (42) 


Equations  (39H42)  are  the  standard  balance  equations  for  the 
number  flux,  momentum,  total  energy  flux,  and  perpendicular 
energy  flux,  respectively. 


As  in  Section  11.  lormal  expressions  lor  may  be  found 
from  (40)  and  (4 1 )  Neglecting  gravity,  we  have 

£, |  =  (2/. *.){«-'  dB/ds(w„i />  -  ..-,.1,) 

+  n"1  d/ds  }  (43) 

£||  =  •«,.)" +  (D/q„nnu„)  d/ds 

■  |(k»Q«||  +  »*o  Qa±)/D\.  (44) 

These  are  the  self-consistent  formulas  for  E h  which  must  be 
satisfied  once  the  solution  for  the  problem  is  found. 

In  order  to  solve  (39H42)  and  (36),  we  must  introduce  clo¬ 
sure  approximations.  Instead  of  using  the  standard  procedure, 
we  introduce  the  same  closure  assumptions  which  we  did  in  Sec¬ 
tion  II-B  by  imposing  (1 7)  and  (18).  Using  (17),  (18),  and  (39), 


we  see  that  the  multimoment  fluid  equations  become 

d/ds(nQua/ B)  =  0  (45) 

d/ds  n0T/;aj|  -  B"1  dB/ds  na[™al(  -  wa±/ 2] 

+  na(qQ/ 2)  d/ds  $ft||  =  0  (46) 

d/ds(wa  j|  +  wa±  +  qQ$  a||)  =  ™0±/uQ  (47) 

B  d/ds{wax/B)  =  was./ua  (48) 

nc  =  ^T  na,  0  =  ions.  (49) 
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Equations  (45H49)  are  a  set  of  4a  +  1  nonlinear,  first-order 
ordinary  differential  equations  for  4a  +  1  unknowns  subject  to 
boundary  conditions  imposed  at  either  end  of  the  flux  tube.  With 
the  aid  of  (45),  we  may  show  that  as  Tq|  0,  (47)  and  (48)  for 
the  ions  reduce  to  the  ion  mean  particle  equations  of  motion 
given  in  [21]  and  [41]. 

The  same  remarks  about  this  closure  assumption  made  in 
Section  II-B  apply  here.  Basically,  this  closure  amounts  to  the 
assumption  that  T0jt  and  mau 2  change  proportionately  from 
si  to  s2,  so  that  the  ratio  (Tajj/2  +  mQu2/2)/(mQu2/2)  is 
nearly  constant.  We  have  examined  some  S3-3  [21],  FREJA 
[22],  and  FAST  [23]  ion  and  electron  particle  data  for  down¬ 
ward  auroral -current  regions.  We  find  that  the  ion  conic  data 
at  the  satellite  altitude  show  distributions  where  both  Tq j  and 
m,ti2  have  increased  more  or  less  proportionately  from  their 
ionospheric  boundary  values.  The  satellite  data  also  show  “flat 
top”  or  “cigar  shaped”  electron  distribution  functions  for  down¬ 
ward  auroral-current  regions  where  a  similar  statement  about 
Te ||  and  meu 2  is  true.  For  these  reasons,  we  argue  that  the  clo¬ 
sure  we  use  is  approximately  valid  for  downward  auroral-cur¬ 
rent  regions. 

Before  we  are  able  to  solve  the  multimoment  fluid  equa¬ 
tions,  we  must  specify  the  turbulent,  resonant  perpendicular  ion 
heating  rate,  ibi±,  and  the  boundary  conditions.  In  this  section, 
we  impose  conditions  consistent  with  the  FREJA  experimental 
observations  by  Boehm  et  ai  [22]  for  a  downward  auroral-cur¬ 
rent  region.  We  do  this  so  that  we  may  compare  the  theoretical 
results  to  the  experimental  results  in  Section  1II-E.  In  that  data 
set,  the  total  electric-field  spectral  density,  S(/),  was  measured 
at  the  satellite  altitude  [22].  In  the  vicinity  of  the  oxygen  cy¬ 
clotron  frequency,  we  found  that  S(f)  ~  f~° .  Since  S  is  eval¬ 
uated  at  /  =  fi ,  this  leads  to 

SL(f  =  fi)  =  ruS0(»/^)3°  (50) 
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Fig.  5.  The  turbulent  (anomalous),  resonant  perpendicular  ion  heating  rate  for  rjL  =  1/2. 1  /4,  and  1/8  from  (37). 


where 

s0  the  radius  of  the  earth,  is  approximated  by  6400  km; 

r/i  is  the  fraction  of  the  total  spectral  power  in  the  LHCP 
ion  cyclotron  component;  and 

a  is  the  spectral  index. 

For  this  data  set,  we  found  that  50  =  2x  10“5  (V/m)2(Hz)“1 
and  a  =  1.3.  Since  no  information  was  available  about  the  alti¬ 
tude  dependence  of  So,  we  assumed  it  to  be  uniform.  In  Fig.  5, 
we  show  a  graph  of  versus  s  for  =  1/2, 1/4,  and  1/8. 
We  see  that  as  the  ions  are  heated  and  move  up  the  B-field 
line,  the  heating  rate  increases  because  the  ions  are  in  reso¬ 
nance  with  lower  and  lower  frequencies.  We  also  assume  that 
the  dominant  ion  was  atomic  oxygen  and  that  si  —  s0  =  800  km, 
s2  -  s0  =  1700  km  (the  satellite  altitude),  £(s)  =  B0($0/&)3, 
and  B0  =  0.57  Gauss.  The  boundary  conditions  at  si  are  n*i  = 
ne i  =  m  =  103  cm’3,  u,i  =  0.5  km/s,  ut i  =  15  km/s, 
Ti in  =  TiJLl  =  0.1  eV,  Te||i  =  Te±l  =  0.3  eV,  and  <t>u  =  0. 
These  parameters  give  jp  =  -2.3  /iA/m2,  w^i  =  0.071  eV, 
and  tuep  =  0.151  eV,  and  are  typical  for  the  topside,  high-Iati- 
tude  ionosphere  [42]. 

We  now  proceed  to  solve  (45H49)  by  a  numerical  method. 
The  method  finds  the  4a  H- 1  values  of  not  ua,  iua|),  iua_Lt  and 
<t>W  as  a  function  of  s  given  their  values  at  one  boundary.  In  this 
case,  si  <  s  <  S2,  and  the  boundary  is  si.  The  procedure  im¬ 
plements  the  semi-implicit  discretization  method  of  Bader  and 
Deuflhard  [43]  and  works  for  nonlinear  systems  of  first-order 
differential  equations  that  may  be  stiff.  The  method  was  tested 
by  setting  Wi±  very  small  and  thereby  obtaining  the  analytical 
results  given  in  Section  II- B  for  those  boundary  conditions.  A 
graph  of  the  numerical  solution  is  given  in  Fig.  6,  where  we  have 
chosen  t)l  =  1/2  and  where  s  measures  the  distance  along  B 
from  the  surface  of  the  earth.  Remember  that  tuax  —  TQ±.  We 


have  also  plotted  the  current  density  j\\ .  Here,  we  have  used  the 
superscript  I  to  denote  the  fact  that  the  particles  are  of  iono¬ 
spheric  origin.  The  dashed  curve  in  the  graph  gives  the 
density  of  the  background  ionosphere  when  there  is  no  turbu¬ 
lent  perpendicular  ion  heating,  i.e.,  when  W{±  =  0-  We  note 
here  that  there  is  a  predicted  potential  increase  of  610  V  from 
800  to  1700  km  for  the  Boehm  et  al.  data  set. 

In  Sections  III-C-III-E,  we  will  discuss  further  aspects  of  the 
solution  shown  in  Fig.  6.  However,  we  wish  to  say  here  that  the 
major  point  of  this  paper  is  that  for  downward  auroral-current 
regions,  E\\  is  negative  and  there  is  a  potential  increase  from 
the  ionosphere  to  the  magnetosphere.  This  is  opposite  to  the  re¬ 
sult  for  upward  auroral -current  regions,  where  the  A-F  formula 
gives  E\\  positive  and  a  potential  drop  from  the  ionosphere  to 
the  magnetosphere.  We  also  note  here  that  is  concave  up¬ 
ward  (dVu/ds2  >  0)  and  that  the  polarization  charge  density 
found  from  Poisson’s  equation  has  the  correct  sign  for  a  nega¬ 
tive  E\\  region. 

C.  Analogue  for  the  A-F  Formula  for  £||  in  Downward 
Auroral-Current  Regions 

Equations  (43)  and  (44)  give  four  self-consistent  formulas  for 
£j |  in  an  electron-ion  plasma.  If  we  had  an  analytic  solution  of 
the  problem,  we  could  substitute  it  into  any  one  of  these  ex¬ 
pressions  and  get  the  correct  formula  for  E\\.  Since  we  do  not, 
we  may  substitute  the  numerical  solution  into  these  equations 
in  order  to  understand  how  a  negative  E\\  is  sustained  in  down¬ 
ward  auroral-current  regions. 

For  example,  consider  (43)  for  the  ions.  It  is 

E\\  =  (2/e){£_1  dB/ds(wi±/2  -  iytj|) 

+  nf 1  d/ds  n,*ttf£|| }.  (5 ! ) 
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Fig.  6.  The  numerical  solution  for  the  turbulent  (anomalous)  heating  of  ionospheric  ions  in  a  downward  auroral -current  region  observed  by  the  FREJA  satellite 
from  (45H49).  The  symbols  are  defined  in  the  text,  rjL  =  1/2,  and  the  superscript  /  denotes  particles  of  ionospheric  origin. 


We  may  integrate  (48)  from  5 x  to  s  to  obtain  a  forma!  expression 
for  Wi± 


Wij_(s)  =  w«x  1  +  B(s) 


ds' 


B(s')m(s')  * 
(52) 


Now  we  can  see  what  is  happening  in  (51).  If  the  turbulent  per¬ 
pendicular  ion-heating  rate  is  sufficiently  large,  then  Wi±/2  - 
ti/j|l  is  positive  (see  Fig.  6)  and,  since  B“l  dB/ds  is  negative, 
the  first  term  on  the  right-hand  side  of  (5 1 )  is  negative.  From 
Fig.  6,  we  also  see  that  n~ 1  d/ds  n,'ii>;||  is  positive,  but  smaller 
than  the  magnitude  of  the  first  term  of  (5 1 )  inside  the  brackets. 
Therefore,  E|j  is  negative.  This  is  shown  in  Fig.  7,  where  Ejjl) 


and  Ejj2)  denote  the  first  and  second  terms  on  the  right-hand 

side  of  (51),  respectively,  and  Ejjtotal)  is  the  total  E\\.  We  see 
that  once  the  numerical  solution  is  known,  then  the  fact  that  E\\ 
is  negative  can  be  understood  from  the  ion  moments  alone. 

For  comparative  purposes,  we  will  refer  to  (51)  and  (52)  as 
the  analogue  of  the  A-F  formula  for  E||  in  downward  auroral- 
current  regions.  Although  (51)  is  formally  the  same  as  (33)  in 
Section  II-C,  its  physical  content  is  different.  For  downward 
auroral-current  regions,  it  is  the  turbulent,  resonant  perpendic¬ 
ular  ion  heating  that  produces  the  ion  temperature  anisotropy 
(wi±  /2  >  )  which  is  consistent  with  the  negative  E\\ .  Since 

upward  and  downward  auroral -current  regions  have  opposite 
ion  velocity-space  anisotropies,  their  parallel  E-fields  have  op¬ 
posite  signs. 
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Fig.  7.  The  solution  of  the  analogue  of  the  A-F  formula  for  £„  in  a  downward  auroral-current  region  given  by  (5 1 ).  Here.  Ej*  ’  is  the  negative  contribution.  EJ” 
is  the  positive  contribution,  and  is  the  total. 


D.  Understanding  the  Behavior  of  the  Solution 

In  Section  III-C,  we  gave  what  we  call  an  analogue  of  the 
A-F  formula  for  E\\  in  downward  auroral -current  regions  and  a 
condition  on  the  moments  of  the  ion  distribution  function  which 
is  self-consistent  with  the  negative  E\\.  Now,  we  wish  to  give  a 
physical  interpretation  of  the  behavior  of  the  ion  and  electron 
moments  shown  in  Fig.  6,  using  the  conservation  relations. 

In  order  to  understand  the  behavior  of  Tjx  and  Tex,  we  recall 
that  wQ±  =  Tax  and  use  perpendicular  energy  balance  (48). 
For  the  ions,  we  have 

d/ds{TiJB)  =  WiJmB .  (53) 

Since  wi±/uiB  is  positive,  Ti±/B  increases  with  s  in  such  a 
way  that  Ti±  increases  from  its  boundary  value  of  0. 1  eV  at  800 
km  to  3.4  keV  at  1700  km.  If  other  perpendicular  ion  heating 
processes  operate  on  the  flux  tube  due,  for  example,  to  the  pres¬ 
ence  of  fast  bipolar  solitary  structures,  then  the  heating  rate  for 
that  process  may  be  added  to  Wi±.  The  amount  of  perpendicular 
ion  heating  would  then  be  increased  from  that  produced  by  the 
BBELF  turbulence  alone.  For  the  electrons,  ti/cx  is  zero,  so 

d/ds(TcjL/B)  =  0  (54) 

and  Tcx  cools  adiabatically  from  its  boundary  value  of  0.3  eV 
at  800  km  to  0.21  eV  at  1700  km.  If  other  electron  heating  pro¬ 
cesses  operate  on  the  flux  tube,  this  result  could  be  modified. 

We  saw  in  Section  III-C  that  when  tUix  is  sufficiently  large,  it 
dominates  the  particle  dynamics  and  the  gravitational  potential 
can  be  neglected  compared  to  the  electrostatic  potential.  Under 
these  conditions,  consider  total  energy  balance  (47)  for  the  elec¬ 


trons.  Since  d/ds  wc±  is  small  compared  to  d/ds  we y,  we  may 
neglect  it  to  find  that 

d/ds  we,|  =  d/ds[Te ||/2  +  meu\t 2]  =  -e£,| .  (55) 

Since  E\\  is  negative,  ujc||  increases  with  s  from  its  boundary 
value  of  0.15  eV  at  800  km  to  610  eV  at  1700  km.  We  see 
that  it  is  the  downward-pointing  £j|  that  is  self-consistent  with 
the  upward  acceleration  of  the  electrons.  Because  of  the  closure 
assumption,  which  is  consistent  with  the  electron  experimental 
data,  we\\  is  partitioned  proportionately  between  the  thermal  part 
and  the  drift  part  so  that  both  increase  with  s.  Equation  (55)  may 
be  integrated  from  si  to  s  to  obtain  an  approximate  conserva¬ 
tion  relation  for  <j>\\  and  tuej|. 

e[<t>^(s)  -  ^j|(si)]  =  u/c||(s)  -  u;c{|(fli).  (56) 

It  is  interesting  to  consider  the  more  rigorous  conservation  of 
total  energy  flux  for  the  electrons  from  (41)  before  any  closure 
approximation  is  made.  Since  d/ds  neQcx  is  small  compared 
to  d/ds  neQc\\  and  the  gravitational  potential  may  be  neglected 
compared  to  the  electrostatic  potential,  we  obtain 

d/ds{(neuJB)[neQcl\/neur  -  e^,]}  ^  0.  (57) 

Using  (39),  we  may  simplify  this  expression  and  integrate  from 
si  to  s  to  obtain 

e[(j>i\(s)  -  <£j|(si)]  =  £o(«)  *“  £o(si)  (58) 

where  E0  is  the  characteristic  energy  and  is  defined  as 
ncQe\\/ncuc.  This  is  just  the  conservation  relation  calculated 
from  the  experimental  data  for  downward  auroral-current 
regions  discussed  and  illustrated  in  Fig.  2  of  Carlson  et  ai 
[23].  When  we  compare  (56)  to  (58),  we  see  that  our  closure 
approximation  is  valid  for  electron  distributions  that  are  more 
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beam-like  lhan  ihermal-like  and  l  here  lore  underestimates  the 
value  ol>u  compared  to  the  more  rigorous  result  given  by  (58). 

Consider  total  energy  balance  (47)  for  the  ions,  where  the 
gravitational  potential  may  be  neglected  compared  to  the  elec¬ 
trostatic  potential 

d/ds[Ti\\/2  +  m{uf/2]  +  d/ds  Ti±  +  r;  d/ds  4>\\  ~ 

(59) 

From  Figs.  5  and  6,  we  find  that  ui{±/ui  is  larger  than 
t  d/ds  <j> || ,  which  is  positive.  Therefore,  the  perpendicular 
wave-ion  heating  rate  supplies  the  energy  to  heat  the  ions 
perpendicular  to  B  and  causes  the  ions  to  be  heated  and  accel¬ 
erated  parallel  to  B ,  as  shown  in  (59).  Because  of  the  closure 
assumption,  which  is  consistent  with  both  the  ion  conic  data 
and  the  kinetic  solution  for  the  ion  “pressure  cooker”  effect 
[29],  i4i||  is  partitioned  proportionately  between  the  thermal 
part  and  the  drift  part  so  that  both  increase  with  s. 

The  rapid  decrease  of  n  can  be  understood  from  conservation 
of  number  flux  (45) 

d/ds(nua/B)  =  0.  (60) 

Since  if,  and  uc  both  increase  with  s  and  B  decreases  with  .9, 
then  n  decreases  rapidly  from  103  cm”3  at  800  km  to  1 1  cm"3 
at  1700  km.  This  produces  a  cavity  in  the  downward  auroral- 
current  region.  For  comparative  purposes,  the  dashed  curve  in 
Fig.  6  gives  the  background  value  of  n  when  there  is  no  turbulent 
ion  heating. 

The  current  density 

jjl  =  en(ui  -  ue)  =  - enue  (61) 

is  dominated  by  the  electron  drift  and  is  therefore  downward¬ 
pointing. 

E.  Comparison  to  a  FREJA  Satellite  Observation  of  a 
Downward  Auroral-Current  Region 

In  this  section,  we  compare  the  theoretical  results  given  in 
Section  IH-B  to  the  experimental  results  of  Boehm  et  al.  [22] 
for  a  FREJA  satellite  pass  through  a  downward  auroral-current 
region  near  1700  km  in  the  morning  auroral  zone.  In  the  Boehm 
et  al.  data  set,  the  total  electric-field  spectral  density,  S(/),  was 
determined  experimentally  at  the  satellite  altitude  and  is  given 
in  Section  III-B.  In  the  theory  we  present,  there  is  one  adjustable 
parameter,  r\i ,  which  we  have  chosen  as  1/2.  Also,  we  assume 
typical  values  [42]  for  the  electron  and  ion  moments  at  the  lower 
boundary,  which  we  chose  to  be  800  km.  A  comparison  between 
the  theoretical  and  experimental  values  for  the  quantities  deter¬ 
mined  by  the  satellite  at  1700  km  is  shown  in  Table  1.  The  ex¬ 
perimental  values  are  taken  from  [22]  where  tic  was  estimated 
using  the  formula  &  - enue ,  and  0  is  the  ion  conic  angle. 
The  other  theoretical  quantities  at  1700  km  may  be  determined 
from  Fig.  6.  We  see  reasonable  agreement  for  all  values  given 
in  the  table  with  the  possible  exception  of  7*  j_.  However,  the  es¬ 
timated  experimental  value  of  5  keV  is  probably  too  large,  and 
a  more  realistic  range  is  from  2  to  3  keV.  It  is  important  to  note 
that  using  the  experimental  value  for  uc  of  940  km/s,  we  find 
that  the  drift  part  of  wc\\  is  2.5  eV  at  1700  km.  From  data  shown 
in  [22,  Fig.  1  ],  we  also  find  that  Tej|  /2  >  meu* /2  at  the  satellite 
altitude.  Apparently,  the  electrons  do  not  run  away  in  an  E\\  field 


.as  a  beam,  but  appear  as  a  field-aligned  distribution  with  a  large 
7'r||/2  compared  to  itt.,  nf  /2  This  is  also  what  is  observed  by 
the  FAST  satellite  |23).  As  pointed  out  in  |22|.  this  is  probably 
due  to  strong  wave-electron  interactions,  it  is  important  to  note 
here  that  wave-electron  interactions  are  not  taken  into  account 
in  the  present  version  of  our  theory,  since  we  use  the  Vlasov 
equation  for  the  electron  dynamics.  Instead,  the  experimental 
fact  that  Tc|j/2  >  m€uf  / 2  is  modeled  by  the  electron  closure 
approximation  that  we  use. 

IV.  Discussion 

In  this  paper,  we  give  an  alternative  derivation  of  the  A-F  for¬ 
mula  for  a  positive  E\\  in  upward  auroral-current  regions  [see 
(30)],  and  propose  an  analogous  formula  for  a  negative  Ejj  in 
downward  auroral -current  regions  [see  (5 1 )  and  (52)].  In  addi¬ 
tion,  we  derive  a  third  result,  i.e.,  an  equivalent  expression  for 
the  A-F  formula  for  £|j  in  terms  of  the  ion  moments  alone  [see 
(33)]. 

The  A-F  formula  is  a  benchmark  result  not  only  for  peda¬ 
gogical  reasons  (it  is  a  closed-form  solution),  but  also  because 
it  clarifies  the  essential  physics  of  a  particular  model  for  au¬ 
roral  arc  formation.  In  that  model,  the  driver  is  the  earth  ward¬ 
streaming,  energetic  magnetospheric  particles,  and  the  mecha¬ 
nism  for  sustaining  a  positive  E\\,  thereby  accelerating  the  au¬ 
roral  particles,  is  the  differential  velocity-space  anisotropy  be¬ 
tween  the  ions  and  the  electrons  in  the  dipolar  magnetic  field 
(tUj||2tt>cX2  -  We|j2^t±2  >  0)-  Another  way  of  looking  at  this 
result,  in  terms  of  the  ion  moments  alone,  reveals  that  E\\  is  pos¬ 
itive  when  wqj  >  Wi±/2  and  certain  other  conditions  are  met. 
According  to  the  PIC  simulation  by  Schriver,  when  ionospheric 
particles  and  wave-particle  interactions  in  the  electrostatic  ap¬ 
proximation  are  included  in  the  model,  the  potential  drop  is  re¬ 
duced  from  its  A-F  value,  and  the  shape  of  <f>\\  is  modified.  In 
his  review,  Borovsky  [17]  identifies  ten  drivers  and  12  possible 
mechanisms,  none  of  which  appears  to  explain  all  aspects  of 
aurora]  arc  formation.  It  is  clear  that  more  analytical  and  simu¬ 
lations!  work  on  upward  auroral-current  regions  is  needed. 

The  primary  result  that  we  present  in  this  paper  is  a  model 
for  a  negative  E\\  in  downward  auroral -current  regions.  The 
BBELF  turbulence,  which  experiment  shows  is  present,  plays 
a  central  role.  The  driver  of  the  turbulence  is  unspecified,  but 
the  spectrum  of  the  turbulence  is  given  experimentally.  The  self- 
consistent  E\\  is  sustained  by  a  two-fold  mechanism  where  the 
BBELF  turbulence  resonantly  heats  the  ions  transverse  to  B 
which,  in  a  dipolar  magnetic  field,  produces  an  ion  velocity- 
space  anisotropy  where  Wi±/2  >  [see  (51)  and  (52)]. 
This  is  just  the  reverse  anisotropy  from  that  which  occurs  in  the 
equivalent  A-F  formula  for  the  ions  [see  (33)]  and  is  consis¬ 
tent  with  an  E y  of  the  opposite  sign.  Of  course,  E\\  and  hence 
<f> j|  must  also  be  self-consistent  with  the  upward  acceleration  of 
the  electrons  [see  (56)  and  the  more  rigorous  expression  given 
by  (58)].  In  this  paper,  we  do  not  attempt  to  resolve  how  the 
downward  auroral-current  region  is  driven.  This  subject  is  an  ac¬ 
tive  area  of  research  and  beyond  the  scope  of  the  present  work. 
For  example,  the  BBELF  could  be  driven  by  the  electron  beam, 
shear  in  the  two-dimensional  ion  flow,  or  some  other  mecha¬ 
nism.  Whatever  the  case  may  be,  the  conservation  laws  derived 
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TABLE  I 

Ihiory-Daia  Comparison  h>k  a  Downward  At  rorai  -Curri  m  Ki.oion 


Quantity  at  1700  km 

Theory 

FREJA  Data 

n  (cm'J) 

11 

-10 

ut  (km/s) 

955 

-940 

ju  (nA/mJ) 

-1.6 

-  -1.5 

n.feV) 

0.21 

0.1  -0.25 

M keV) 

3.4 

<5 

Mb  (Volts) 

610 

-700 

6  (degrees) 

106 

. . 

100-110 

from  the  appropriate  kinetic  model  impose  powerful  constraints 
on  the  moments  of  the  electron  and  ion  distribution  functions, 
the  strength  of  the  BBELF  turbulence,  and  the  magnitude  and 
direction  of  E\\.  These  conservation  relations  were  discussed  in 
Sections  III-B  and  III-D.  The  central  point  of  this  paper  is  that 
both  the  velocity-space  anisotropy  of  the  particles  and  the  tur¬ 
bulent  heating  of  the  ions  play  an  important  role  in  downward 
auroral-current  regions. 

When  this  model  for  downward  auroral -current  regions  is 
applied  to  the  FREJA  data  set  by  Boehm  et  ai  [22]  in  the 
morning  auroral  zone,  good  agreement  between  theory  and 
experiment  is  achieved  (see  Table  I).  For  example,  we  find  that 
jy  ~  “1.6  jzA/m2  at  1700  km,  the  average  parallel  £-field 
<£|j)  ~  —680  /zV/m,  and  a  potential  increase  of  ~610  V 
from  800  to  1700  km  are  obtained.  In  [29]  and  {30],  a  kinetic 
theory  of  the  “pressure  cooker”  effect  was  given  for  the  S3-3 
data  set  of  Gomey  et  ai  [21]  using  (34M36)  for  a  downward 
auroral -current  region  in  the  dayside  auroral  zone.  Again,  good 
agreement  between  theory  and  experiment  was  achieved.  Since 
the  turbulence  was  weaker  than  that  in  the  FREJA  data,  we 
found  that  j?||  ~  —0.4  /zA/m2  at  6000  km,  {JE||)  ~  — 12  /xV/m, 
and  a  potential  increase  of  ~55  V  from  1500  to  6000  km 
were  obtained.  In  the  anomalous  transport  calculation  for  the 
current-driven  ion  cyclotron  instability  reported  in  [26],  they 
found  that  j\\  ^  — 1.1  /zA/m2  and  (J£||)  ~  -0.1  /zV/m  for  a  1 
V  potential  increase  in  104  km.  This  result  for  (£j|)  appears  to 
be  loo  small  to  be  of  significance  in  downward  auroral-current 
regions. 

The  model  presented  here  for  downward  auroral -current  re¬ 
gions  needs  to  be  generalized.  There  is  the  obvious  need  to  iden¬ 
tify  the  driver  of  the  BBELF  turbulence  and  to  treat  the  wave- 
particle  interactions  self-consistently.  Aside  from  this,  there  are 
other,  more  straightforward  generalizations  that  can  be  carried 
out.  Recently,  we  have  generalized  this  work  in  two  ways.  In 
addition  to  the  presence  of  the  heated  and  accelerated  iono¬ 
spheric  particles,  we  have  included  a  magnetospheric  popula¬ 
tion  of  earth  ward- streaming  ions  and  electrons.  We  find  that  as 
the  flux  of  the  magnetospheric  particles  is  increased,  a  positive 
contribution  to  E\\  is  produced,  and  the  total  E\\  is  less  negative, 
in  fact,  if  the  flux  of  magnetospheric  particles  is  large  enough, 
the  two  effects  combine  to  produce  a  positive  E n  and  an  upward 


auroral-current  situation.  A  second  way  in  which  wc  have  gener¬ 
alized  the  model  is  to  include  wave-particle  interactions  in  the 
equation  for  the  electrons.  The  BBELF  waves  cannot  heat  the 
electrons  perpendicular  to  B  hy  a  cyclotron  resonant  interaction 
because  the  frequencies  are  too  low.  but  there  is  a  nonresonant 
diffusion  process  which  acts  parallel  to  B.  Preliminary  results 
indicate  that  this  diffusion  process  spreads  out  the  electron  dis¬ 
tribution  function  in  velocity  space  parallel  to  B ,  prevents  the 
formation  of  a  runaway  electron  beam,  and  maintains  a  parallel 
electron  energy,  t/;c||,  where  TtfJ|/2  >  mr?z2/ 2. 

It  should  be  mentioned  here  that  there  are  processes  other 
than  ion  cyclotron  resonance  with  broadband  waves  by  which 
ions  can  be  heated  in  space  plasmas.  For  a  discussion  of  these 
processes,  see  [44).  Recently,  a  new  process  has  been  proposed 
based  on  an  analysis  of  FAST  data.  The  data  show  that  within  a 
broad  region  of  downward  auroral  current,  there  are  narrow  sub- 
regions  where  intense  bursts  of  fast  bipolar  solitary  structures 
occur.  These  subregions  coincide  with  enhanced  ion  heating 
(see  [24,  Fig.  1]).  It  has  been  pointed  out  [23]  that  the  random 
pulses  in  the  perpendicular  £-field  of  the  solitary  structures  as 
they  propagate  past  the  nearly  stationary  ions,  could  lead  to  a 
significant  amount  of  stochastic  ion  heating.  If  a  perpendicular 
ion  heating  rate  can  be  calculated  for  this  process,  then  it  may  be 
added  to  Wi±  in  (41)  and  (42).  The  resulting  multimoment  fluid 
equations  may  then  be  solved  using  the  same  method  as  pre¬ 
sented  in  Section  III-B.  This  would  lead  to  an  enhanced  amount 
of  perpendicular  ion  heating,  as  indicated  in  [24,  Fig.  1]. 
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